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Experiments performed at high gas partial pressures have demonstrated that the kinetics of the 



CO oxidation reaction at Ru (0001) is different and somewhat anomalous compared to that over 
other transition metal surfaces and, in particular, the turnover rate is exceptionally high. In or- 



J> ■ der to gain insight into the underlying reasons for this behavior, we performed density functional 



theory calculations using the generalized gradient approximation for the exchange-correlation 
functional. We find that the high rate is due to a weakly, but nevertheless well bound, (1 x 1) 
oxygen adsorbate layer which may form for high O2 pressures but not under usual ultra high 



vacuum conditions. The calculations indicate that reaction to CO2 occurs both via scattering of 



■ gas-phase CO molecules as well as by CO molecules weakly adsorbed at vacancies in the oxygen 

C . 

O ■ adlayer, where the latter mechanism dominates the rate. 

> 

53 . I. INTRODUCTION 

The oxidation of carbon monoxide at transition metal surfaces is one of the most ex- 
tensively studied and best understood heterogeneous catalytic reactions (see, for example 
Refs. jlj and references therein). Nevertheless, very little in fact, is actually known about 
the reaction pathway on a microscopic level. Recently there has been a resurge of interest 
in such surface chemical reactions which has been stimulated by advances in surface science 
techniques (e.g. molecular and atomic beam experiments and high pressure catalytic 
reactors ||). These techniques have enabled new information to be obtained concerning the 
behavior of chemisorbed reaction partners. Recent studies using high gas pressure catalytic 



reactors (e.g. operating at ~ 10 torr with CO/O2 pressure ratios < 1) have reported unusual 
behavior for the CO oxidation reaction over Ru(0001) Q: (i) The rate of CO2 production 
is significantly higher than at other transition metal surfaces; in contrast, under ultra high 
vacuum (UHV) conditions, Ru(0001) is notably the poorest catalyst for this reaction [|TJ. 
(ii) The measured kinetic data (e.g., activation energy, and temperature and pressure de- 
pendencies of the rate) is markedly different to that of other substrates. (Hi) As opposed to 
other transition metal catalysts (Pt, Pd, Ir, and Rh), highest rates of CO2 formation occur 
for high surface oxygen concentrations (one monolayer was proposed on the basis of auger 
electron spectroscopy measurements), (iv) Almost no chemisorbed CO was detected during 
or after the reaction. On the basis of these findings, it was speculated that an Eley-Rideal 
(E-R) mechanism is operational as opposed to the "usual" Langmuir-Hinshelwood (L-H) 
mechanism. In the E-R mechanism, the reaction occurs between gas-phase and chemisorbed 
particles, while in the L-H mechanism, the reaction occurs between species both chemisorbed 
on and in thermal equilibrium with the surface. So far E-R mechanisms have only recently 
been experimentally confirmed, and in most cases they are for somewhat artificial reactions 
between gas-phase atomic hydrogen or deuterium (from a beam source) with chemisorbed 
atoms jS|. 

In the present paper we describe results our theoretical study aimed at understanding 
the experimental findings described above for the oxidation of CO at Ru (0001). 

II. CALCULATION METHOD 

We performed density functional theory (DFT) calculations using the pseudopotential 
plane wave method and the supercell approach 0. Surfaces were modelled by unit cells of 
(2 x 2) and (1 x 1) periodicity with four layers of Ru (0001) and a vacuum region equivalent 
to thirteen such layers. The ab initio fully separable pseudopotentials were created using the 
scheme of Troullier and Martins 0. Plane wave cut-offs of 40 Ry, with three k-points || in 
the surface Brillouin zone of the (2 x 2) cell were used, as well as a 60 Ry cut-off with fifty-two 
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k-points in the surface Brillouin zone of the (lxl) cell. Results of calculations performed 
using these different basis sets indicated that the smaller basis set provides a sufficiently 
accurate description for the present investigation ||. The generalized gradient approxima- 
tion (GGA) of Perdew et al. [IIJ was employed for the exchange-correlation functional and 
it was used in creating the pseudopotential as well as in the total-energy functional. It is 
therefore treated in a consistent way. The position of all atoms were allowed to relax except 
the bottom two Ru layers which were held at their bulk positions. We like to point out 
that the present study represents the first investigation of a surface chemical reaction (to 
be distinguished from dissociative chemisorption) using such theoretical methods. Further 
details about the calculations can be found in Refs. |P, |TT| , |T2" . 



III. ADSORPTION OF O ON RU (0001) 

Given that the CO oxidation experiments identified that highest rates of CO2 produc- 
tion occurred for high O2 partial pressures and concomitantly high oxygen coverages, an 
important part of the theoretical study is therefore to investigate the structure and stability 
of high coverage O adlayers on Ru (0001). It is known that using O2, under UHV conditions 
at 300 K, the saturation coverage is close to 8 = 0.5. Earlier we reported the results of 
DFT-GGA calculations for various adsorbate phases of oxygen on Ru (0001) and found that 
an even higher coverage should be stable on the surface ||. In particular, a structure with 
a (1 x 1) periodicity and coverage = 1 where the oxygen atoms occupy the hep-hollow 
sites. (The coverage O is defined as the ratio of the concentration of adparticles to that of 
substrate atoms in the topmost layer.) This result indicated that the formation of the = 1 
structure under UHV conditions from gas phase O2 is kinetically hindered, but that by of- 
fering atomic oxygen, this phase should be attainable. Indeed this structure was successfully 
created by starting from the = 1/2 phase and offering additional atomic oxygen to the 
surface via the dissociation of NO2 ||12|| . The determined surface atomic structure of the 



[1 x 1) phase by the low-energy electron diffraction (LEED) intensity analysis of Ref. [12" 
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agreed very well with that predicted by the DFT-GGA calculations. The calculations show 
that the adsorption energy of oxygen in this phase is significantly weaker than it is in the 
lower coverage ordered phases (of (2 x 2) and (2 x 1) periodicities). This can be seen from 
Fig. 1 which shows the adsorption energy, with respect to O2 gas, for oxygen in the hcp- 
and fcc-hollow sites at various coverages. The region over which dissociative chemisorption 
of oxygen is kinetically hindered is indicated as the shaded region. It can be seen that the 
hep-hollow site is favored at all coverages, where the adsorption energy of the (2 x 2) phase 
is the greatest. This suggests that at low coverages island formation will occur with a (2 x 2) 



periodicity which is consistent with experimental observations [0,0. With subsequently 
increasing coverage, the adsorption energy decreases which reflects a repulsive adsorbate- 
adsorbate interaction, and indicates that no island formation will take place in the coverage 
regime G = 0.25 to 1. Furthermore, with increasing coverage from Q = 0.25, the difference 
in adsorption energy between the fee- and hep-hollow sites decreases, and at coverage one 
the difference is very small (~0.06 eV). 

Catalysis involving high O2 pressures is therefore likely to involve oxygen coverages that 
approach one monolayer on the surface because the impingement rate is proportional to 
the partial pressure which means that there is a significantly higher attempt frequency to 
overcome activation barriers for dissociative adsorption than under UHV conditions. We 
therefore initially assumed in our investigation of the oxidation of CO at Ru (0001) that the 
(lxl) phase covers the surface. 



IV. OXIDATION OF CO AT (1 x l)-0/RU (0001) 

For the catalytic oxidation of CO, it is necessary that a CO molecule reacts with an oxy- 
gen atom to yield CO2. From the experiments it had been speculated that CO may achieve 
this from the gas-phase without adsorbing on the surface first, i.e. the E-R mechanism. To 
investigate the possibility of this scenario, calculations were first carried out to determine 
whether CO could adsorb on the (1 x 1) oxygen covered Ru(0001) surface. The sites con- 
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sidered were the on-top and fcc-hollow sites, with respect to the Ru (0001) substrate, and 
a bridge site between two adsorbed O atoms. It was found that at all these sites, CO is 
unstable, thus preventing the L-H mechanism on the perfectly ordered (1 x l)-0/Ru (0001) 
surface. 

To investigate the possibility of reaction via an E-R mechanism, we evaluate an appropri- 
ate cut through the high-dimensional potential energy surface (PES); this cut is defined by 
two variables: the vertical position of the C atom and the vertical position of the O adatom 
below the molecule. In order of ease of analysis, the CO axis is held perpendicular to the 
surface. The resulting PES is presented in Fig. 2 where the coordinate system is shown as 
the inset. The most favorable pathway can quickly be seen: it has an energy barrier of about 
1.6 eV and formation of CO2 is achieved via an upward movement of the O adatom by ~ 0.4 
A towards the CO molecule. Due to the similar masses of O and C, it is likely that the 
impinging CO molecule will impart a significant amount of energy to the O adatom, thus 
stimulating its vibrations and facilitating its motion (indicated by the oscillations in the 
dot-dashed line). This would then bring the system to the transition state of the reaction 
(marked by the asterisk). The newly formed CO2 molecule is strongly repelled from the 
surface and moves into the gas phase region with a large energy gain of 1.95 eV. 

As mentioned, the PES of Fig. 2 corresponds to a constrained situation of the surface-CO 
angle. When this constraint is dropped, i.e., when the tilt angle of the CO axis is allowed 
to relax, we find that the energy barrier is reduced to 1.1 eV, and also that the position 
of the saddle point occurs closer to the surface (by 0.3 A). The geometry of the identified 
transition state at the saddle point is depicted in Fig. 3. The optimum tilt angle with respect 
to the surface normal is 49° which corresponds to a "bond angle" of 131° for the "C02-like" 
complex. It is interesting to note that this geometry is similar to that associated with 



the CO2 ion [[HJ and to that proposed for the "activated complex" for the CO oxidation 
reaction over other transition metal surfaces |16[] in which the reaction proceeds via the L-H 
mechanism. 

In Fig. 4 the energy diagram for the proposed reaction pathway for the E-R mechanism 



is given. Interestingly, the calculations show that there is a physisorption well for CO 
(over all of the surface unit cell), as well as for CO2, above the surface. The calculated 
depths however (~0.04 eV), are likely to be underestimated because the employed exchange- 
correlation functional does not describe the long-range van der Waals type interactions. It 
can be seen that with respect to the free CO and |02 there is an energy gain of 1.20 eV on 
adsorption of oxygen into a vacancy in the (1 x 1) structure, an energy barrier for reaction 
to CO2, and a significant energy gain on formation of CO2 which leaves an oxygen vacancy 
behind on the surface. This predicted mechanism and associated energetics could possibly 
be tested by molecular beam experiments. It is, however, unlikely that this process alone 
can explain the high rate measured in the high pressure catalytic reactor experiments given 
the height of the activation energy barrier, which is high for an E-R mechanism. 

It is conceivable that CO molecules adsorb at vacancies in the O-adlayer, i.e., at sites 
at which an oxygen atom has been removed (e.g. by the above described E-R reaction). 
We find that CO can weakly adsorb in such a vacancy where the adsorption energy is 
0.85 eV- significantly less than on the clean surface which we calculate to be 1.68 eV (the 
experimental value has been determined to be 1.66 eV [I8"|). Adsorption of |02 in such a 



vacancy is more favorable, with an adsorption energy of 1.20 eV. Using these adsorption 
energies of CO and |0 2 , and assuming thermal equilibrium of the CO + 2 gas and a 
mixed CO + O adlayer, the law of mass action indicates that about 0.03 % of the sites 
of the (1 x 1) adlayer will be occupied by CO. (We assumed that the O2 and CO partial 
pressures are equal and the temperature is T = 500 K.) In reality the CO concentration 
will be even higher because we find that CO adsorption into an existing O vacancy can 
proceed basically without hindrance while O adsorption (from O2) is hindered by an energy 
barrier. Therefore the actual percentage of surface sites occupied by CO will be somewhat 
higher. A L-H reaction between the adsorbed CO molecule in a vacancy and a neighboring 
O adatom is expected to be particularly efficient and to give rise to the high rate due to the 
weaker bond strengths of both the adsorbed CO molecule and O adatom, as well as their 
close proximity to one another. In contrast, under UHV conditions, where the high coverage 
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phase is unable to form due to kinetic hindrance for 2 dissociation, Ru (0001) is the poorest 
catalyst for this reaction. The reason for this is likely to be related to the fact that under 
these conditions ruthenium binds oxygen (and carbon monoxide) particularly strongly. 



7 



REFERENCES 

[1] T. Engel and G. Ertl, Adv. Catal. 28 (1979) 1; The chemical physics of solid surfaces 
and heterogeneous catalysis, eds. D. A. King and D. P. Woodruff, Vol. 4, Elsevier (1982). 

[2] K. D. Rendulic and A. Winkler, Surf. Sci. 299/300 (1994) 261. 

[3] D. W. Goodman, Surf. Sci. 299/300 (1994) 837. 

[4] C. H. F. Peden, in Surface Science of Catalysis: In situ probes and reaction kinetics, 
Eds. D. J. Dwyer and F. M. Hoffmann, Am. Chem. Soc, Washington DC (1992). 

[5] E. W. Kuipers, A. Vardi, A. Danon, and A. Amirav, Phys. Rev. Lett. 66 (1991) 116; C. 
T. Rettner, Phys. Rev. Lett. 69 (1992) 383; C. T. Rettner and D. J. Auerbach, Phys. 
Rev. Lett. 74 (1995) 4551; C. T. Rettner, J. Chem. Phys. 101 (1994) 1529. 

[6] R. Stumpf and M. Scheffler, Comp. Phys. Commun. 79 (1994) 447. 

[7] N. Troullier and J. L. Martins, Phys. Rev. B 43 (1991) 1993. 

[8] S. L. Cunningham, Phys. Rev. B 10 (1974) 4988. 

[9] C. Stampfl and M. Scheffler, Phys. Rev. B, 54, 2868 (1996). 

[10] J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, 
and C. Fiolhais, Phys. Rev. B 46 (1992) 6671. 

[11] C. Stampfl, Surf. Rev. Lett., in press. 

[12] C. Stampfl, S. Schwegmann, H. Over, M. Scheffler, and G. Ertl, Phys. Rev. Let., in 
press. 

[13] K. L. Rostov, H. Rauscher, and D. Menzel, Surf. Sci. 278 (1992) 62. 

[14] J. Wintterlin, J. Trost, S. Renisch, R. Schuster, and G. Ertl, to be published. 

[15] J. Pacansky, U. Wahlgren, and P. S. Bagus, J. Chem. Phys. 62 (1975) 2740. 

8 



[16] G. W. Coulston and G. L. Haller, J. Chem. Phys. 95 (1991) 6932, and references therein. 

[17] W. J. Mitchell and W. H. Weinberg, J. Chem. Phys. 104 (1996) 9127. 

[18] H. Pfnur, P. Feulner, H. A. Engelhardt, and D. Menzel, Chem. Phys. Lett. 59 (1978) 
481. 



9 



FIGURES 



> 

w 
>^ 

0) 

a-) 
o 

• 1— I 

» 

& 

o 

< 



1.6 
1.8 
2.0 
2.2 
2.4 
2.6 
2.8 





1 i 1 i 1 -t 


ice site 








• 






' kin. 




hindered _ 




regime 


hep site 

ii 


1,1,1 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Oxygen coverage 



FIG. 1. Adsorption energy of O on Ru(0001) with respect to \02 at various coverages for 
oxygen in the fee- (dashed line) and hep-hollow sites (continuous line). 
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FIG. 2. Potential energy surface (PES) as a function of the positions of the C atom, Zq, and 
of the O adatom, Zq. Positive energies are shown as dashed lines, negative ones as full lines. The 
contour line spacing is 0.6 eV. 
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Transition State 




FIG. 3. Illustration of the transition state identified for the reaction of gas-phase CO with 
adsorbed oxygen when tilting of the CO axis is considered. 
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FIG. 4. Predicted energy diagram for the E-R mechanism of CO oxidation at Ru (0001). Note 
that the depths of the physisorption wells are exaggerated for visibility. 
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This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9611217vl 



